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ABSTRACT 

We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and 
an average cosmological field at ^ ^ 6 — 12, in order to investigate the influence of environment 
and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic 
winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this 
evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable 
wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to 
accelerated evolution of dark matter and baryonic structures, absent in the “normal” region, and to 
shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on 
both environments, the more physically motivated VW model does exhibit this effect. In addition, 

VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW 
and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the 
intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas 
which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in 
galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from 
deficiencies, the VW model seems to be promising in correlating the outflow properties to those of 
host galaxies. Further constraints on the state of the IGM at high- 2 : are needed to separate different 
wind models. 

Subject headings: cosmology: dark ages, reionization, first stars — cosmology: theory — galaxies: 
formation — galaxies: high-redshift — methods: numerical 


1. INTRODUCTION 

With the rapidly increasing number of available multi¬ 
wavelength observations as well as the development of 
more sophisticated theoretical and numerical tools, our 
understanding of galaxy formation and evolution has pro¬ 
gressed considerably in the recent years. While the ma¬ 
jority of large galaxy surveys have focused on low and in¬ 
termediate redshifts, substantial data are now available 
to gain insight into the early stages of galaxy formation 
at 2 : > 6. In particular, observations of high-redshift 
quasars (QSOs) at 2 : > 6 have hinted at the presence of 
rare, overdense regions hosting massive str uctures less 
than one billion years after the Big Bang (Fan et al. 
2001). Deep observations with the wide-area buprime- 
Cam of the Subaru telescope have led to the discovery 
of ^ 3<j overdensities of Lyman Break Galaxies (LBGs) 
on scales of a few Mpc around th e 2 : ^ 6.4 QSO CFHQS 
J2329-0301 (Utsumi et al. 2010). Similar results have 
been obtained using the Large Binocular Telescope to 
study the environment of other high-^ QSOs and also 
concluded that they reside in overdense regions (e.g.. 
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Recent theoretical work has shown 
that these regions do not necessarily evolve into the most 
massive clusters at ^ = 0 and argued that galaxy forma¬ 
tion might proceed differently there than in average den¬ 
sity fields dTrenti fc Stiavell i||2QQ8| |Overzier et al.||2009 
Romano-Diaz et aL||2Ulla|b|). 

Star formation and AG activity both appear to pea k 
around ^ ^ 2 — 3 (e.g. Madau & Dickinson 2014), 
and galactic outflows, or galactic winds, are expected 
to exhibit a similar behavior, when averaged over a 
large enough volume. However, overdense regions, such 
as those hosting high- 2 : QSOs, have been predicted to 
evolve ahead of the average and under dense regions (e.g. 


Barkana fc Loeb|2QQ4{|Romano-Dfaz et al.|2Qlla|b[|2Q14[ 


Yajima et al.||2015). Although their evolution is concur- 

rent with the reionization epoch at 6 < 2 : < 14, they 


can be reionized ahead of the normal regions in the uni¬ 
verse because of the elevated star formation rates there. 
Moreover, these regions should collapse early, and, be¬ 
cause the universe is also substantially denser at these 
redshifts, one expects mass accretion rates to be high. 
The question then is to what degree the associated galac¬ 
tic outflows from these objects affect their evolution as 
well as the evolution of their environment. In this paper, 
we aim to study the effects of these winds on galaxy evo¬ 
lution using high-resolution cosmological simulations of 
overdense regions and compare them with the evolution 
in average-density regions. 

Massive outflows are projected to play an important 
role in the evolution of a large number of astrophysi- 
cal objects — from single and binary stars to supermas- 
sive black holes (SMBHs) and starburst galaxies. On 
galactic scales, winds are expected to affect the evolu- 
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tion of the ir hosts by reg:ulating: the star formation rate 
(SFR , e.^.jScannapieco et aL|2QQ6||Oppenheimer fc Dave 
2006 Schaye et ah '20101), ejecting gas and metals into 

tic medium (IGM, e.^., [Tremonti et al.||2004 


et al. 201^. The contribution of AGN is debatable at 
The 


the parent dark matte r (DM) halos and the interg:alac- 

ni[2or" 


Dave et al 


2008 Kirby et al.|20Q8 [bijacki et al.|2009[ ), and by deter 


mining the physical properties (disk mass and size, chem- 


preseiit. I'he driving of the outflows by the SN and stel¬ 
lar winds occurs when ejecta of individual sources form 
a bubble of hot gas, ^ 10^“^ K, which expands due to 
the strong overpressure down the steepest pressure gra¬ 
dient and enters the ‘blow-out’ stage. Besides the mass, 
energy and momentum injected by such winds, this is 


of DM halos) of these objects (e.g., Dekel & Silkll986| 

1 El- to be placed into the halo and further out into the IGM 

Zant et al.||2001 2004 Romano-Diaz et al.||2008p. 'I'l 

liese (e.g., 

Cen et al.||2005). 


outflows powered by supernovae (bW) can be supple¬ 
mented by a feedback from the accretion processes onto 


called AGN winds (e.g., Blandford & Payne|1982 

Shlos- 

man et al.|| 19851 Fmmerin^ 

^ et al.||1992 

Konigl ^ 

Kartje 

1994 

Murray et al.l 1995 | 

[Arav et al. 

19971 IGecil et al. 

2001 

servai 

Proga 2003p for whi 

tional evidence (e.g.. 

ch there exists numerous ob- 

Emonts et al. 2005 Croston 

et al. 120081 ICano-Dfaz et 

al.||2012t IGombes et al. 

ponr 

Cicone et al. 2014p. A compelling example which has 


can be found in the comparison study of disk evolu¬ 
tion in a c osmological setting, wi t h and without feed - 
back (e.g., Robertson et al. 2004 Schaye et al. |2015), 
where in the absence of a feedback the gas quickly vio- 
lates the Toomre criterion and fragments. Without feed¬ 
back there is an overproduction of metals, especially in 
small galaxies. Even more revealing is the over cooling 
problem which leads directly to the angular momentum 
catastrophe — the gas falls into the subhalos, and hitch¬ 
hikes to the bottom of the poten tial well of the parent 
galaxy (e.g.. Mailer fc Dekel|2002 ), possibly contributing 


to the overgrowth of classical bulges. 

The presence of galactic outflows is supported by nu- 


ies (e 

.g., Osterbrock 

1960 Lynds & Sandage|1963l 

Math- 

ews ^ 

V Baker||19711 

ieckman et al.||1990 

|Humme 

et al. 

I991I 

Heckman 1994| Kunth et al.||1998 

), including the 

2 ; ^ 3 

1 — 4 LBGs (e.g.,|Pettini et al.|2000[jAdelberger et al. 

2003 

Shapley et al.|2003|), gravitationally-lensed galaxies 

at 2 ; ^ 4 —5 (e.g.,|hranx et al.|1997), luminous and ultra- 


etc. About 75-100% of ULIRGs are associated with such 
winds — their strength appears to correlate with the star 
formation rates (SF Rs), but saturates for the most lumi - 
nous ULIRGs (e.g., [Martin 2005 Veilleux et al. 2005). 


Galactic winds are typically characterized by a bicomcal 
sym metry with respect to the underlying galactic disks 


(e.g., Veilleux et al. 1994 Shopbell & Bland-Hawthorn 


1998|). Hot tenuous winds show eviden ce tor associated 
10^ K and n eutral clouds (e.g., Heckman 


cold __ 

Stewart et al. 2000| ) 


_ 1994 

so they represent a multiphase 
interstellar medium (ISM). An important link with the¬ 
oretical predictions has been made by recent observations 
revealing that outflow velocities appear to correlate w ith 
the galaxy stellar mass or its SFR (e.g., Martin||2005 ). 


A number of driving mechanisms have been proposed 
for galactic winds, such as radiation pressure in the UV 
lines and on dust, thermal pressure from the SN and OB 


Castor et al. 


Dekel 


& Silk|19861|Mac Low ^ lVlcCray|1988 

lOstriker McKee 

19881 

Shull & Saken|19951 Costa et al. 

2014 Vogelsberger 


understanding of various mechanisms that regulate the 
star formation and angular momentum redistribution in 
a forming disk galaxy. An unexpectedly large fraction, 
^ 76%, of massive, > 10^^ Mq galactic disks can be fit 
with the Sersic index < 2, and ^ 69% have a bulge-to- 
the-total disk mass (B/T) r atio of 0-2, both in barred 
and unbarred galaxies (e.g., [Weinzirl et al. 2009). Such 


low Sersic indexes repres ent disky rather then spheroida l 


stellar distribution (e.g., [Kormendy fc Kennicirtt 


2004). 


simu- 


On the other hand, bulges obtained in numerica. 
lations are dominated by massive spheroids. This con¬ 
tradiction becomes even stronger in the recent study of 
the bulge population within a sphere of llMpc radius 
around th e Milky Way using Sp itzer 3.6 /im and HST 
data (e.g., Fisher & Drory||2011 ). The dominant galaxy 
type in the local universe has been found to possess pure 
disk properties, i.e., having a disky bulge or being bulge¬ 
less. These results reinforce the opinion that additional 
physical processes are needed to explain much less mas¬ 
sive spheroidal components in disk galaxies. 

Overall, strong arguments exist in favor of a pro¬ 
cess which lowers the effici ency of gas-to-stars conversion 
(e.g., Fukugita et al.|1998| ). Such a process (or processes) 
may resolve the discrepancy between the observed galaxy 
luminosity function (LF), both for high- and low-mass 
galaxies, and the computationally-obtained DM Halo 
Mass Function (HMF). Simulations of baryon evolution 
within the A Cold DM (ACDM) framework appear to 
overproduce low and high-mass galaxies. Feedback from 
stellar evolution in low-mass galaxies and from AGN in 
high-mass ones is expected to quench the star formation, 
depleting the host galaxy from its interste llar medium 
an d resolving this discrepancy (e.g., [Khochfar fc 


Carthy et al.||2012| |Hi Iz et A ■ II2013II , see al so reviews w 

Veilleux et al.l (|2U05j) and Shlosman (2013). 

In this work we focus at rare overdense regions (by 
construction) in the universe which host massive ^ 
10^^ h~^ Mq dm halos by 2 ; ^ 6, possibly the host halos 
of QSOs, and quantify the evolution of DM and baryons 
in these regions, in comparison with the ‘normal’, aver¬ 
age density regions. We compare numerical models of 
galactic winds by means of high-resolution cosmological 
simulations of galaxy formation and evolution at these 
high redshifts, and analyze efficiency of these winds in 
carrying metals outside galaxies, into the halo and IGM 
environment. The effect of these winds on individual 
galaxies is discussed as well. 

This paper is structured as follows. Section 2 describes 
the numerical methods, initial conditions used, and the 
details of the wind models. Resulting halo and galaxy 
properties are given in section 3, and the wind effects 
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on the large-scale environment are analyzed in Section 4. 
Discussion and conclusions are given in the last section. 

2. SIMULATIONS 
2 . 1 . The code 


the method of SHQ3 and a variable ve locity wind (VW) 
model from Ghoi & Nagamine (2011 hereafter GNU), 
supplemented with a model without outflows (NW). All 
of our three wind models share the same thermal feed¬ 
back as in SH03. 


We use a modified version of the tree-particle- 
mesh Smoothed Particle Hydrodyn amics (SPH) c ode 
GADGET-3, originally described in Sprmgel (|2QQ5|), in 


quist 


its conserva tive entropy formulation (Springel Hern- 
' Our conventional code includes radiative 


2002J 


cooling by H, He, and metals (|Ghoi & Nagamine|2009|), a 
recipe for star formation and SN feedback, a phenomeno- 
logical model for galactic winds, and a sub-resolution 


mode l for the multiphase ISM (Springel & Hernquist 


2003 hereafter SH03). In the rnultiphase ISM model. 


star forming SPH particles contain the cold phase that 
forms stars and the hot phase that results from SN heat¬ 
ing. The cold phase contributes to the gas mass, and 
the hot phase contributes to gas pressure. Metal enrich¬ 
ment is also taken i nto account according to the recipe 
from SH03 (see also Choi & Nagamine 2009) in which 
the metallicity increase of star-forming gas particles is 
related to the the fraction of gas in the cold phase, the 
fraction of stars that turn into SN and the metal yield 
per SN explosion. Although metal diffusion is not explic¬ 
itly implemented, metals can still be transported outside 
galaxies by wind particles and enrich the halos and the 
IGM. 

Since we only focus on 2 : > 6 , before the full reioniza¬ 
tion, and do not implement on-the-fly radiative transfer 
of ionizing photons, the UV background is not included 
in our simulations. One expects this omission to have a 
larger effect on the overdense regions. To quantify this 
would require an introduction of additional free param¬ 
eter (s), as the precise level of the UV background is un¬ 
known at present. We, therefore, refrain from doing so. 
Eor the sa me reason we have neglected the AGN feed¬ 
back (e.g., Sijacki et al. 2009) which has been predicted 
to affect the more massive galaxies. 

Eor the star formation, we use the “Pressure model” 


(Choi & Nagamine 2010) which implements the rela- 
tionship between the gas surface density and its vol¬ 
ume density using |Schaye & Dalla Vecchia| (|2008|) pre¬ 
scription. The characteristic time-scale tor the bf' be¬ 
comes rsF ^ A“^(l Mq pc“^)^ whnre 


A = 2.5T 0.7M^ yr-\ n = 1.4 ± 0.15, 7 = 5/3 (|Kenni- 
|cutt||1998[ ), and P is the total effective gas thermal pres¬ 
sure (including the contribution from both cold and hot 
phases). This model reduces the high- 2 : SER relative to 
the standard recipe described in SH03. Star formation 
is triggered when the gas densi ty is aboye the thresh¬ 


old value = 0.6 cm ^ (e.g., Springel 


2005 


Choi & 


Nagamine 2010 Romano-Dfaz et al.| 2011 b ) . T he value ot 


is based on the translation of the threshold surface 
density in the Kennicutt-Schmidt law, SER r\j Sgas> where 
SER is the disk surface density of star formation, Egas 
is the surface density of the neutral gas, and a ^ 1 — 2 , 
depending on the tracers used and on the relevant linear 
scales. 


2 . 2 . Galactic wind models 

We have considered three different galactic outflow 
models: a constant velocity wind (GW) model based on 


No-Wind model (NW )— This model has no kinematic 
feedback, but maintains thermal feedback by the SN. 

Constant Wind model (CW )— The galactic wind has been 
triggered by modifying the behavior of some gas parti¬ 
cles into the ‘wind’ particles. Those were not subject 
to hydro dynamical forces and have experienced the ini¬ 
tial kick from the SN. All CW particles had the same 
constant velocity, Uw = 484kms“^, and the same mass¬ 
loading factor, ydw = Mw/Msf, where is the mass 
loss in the wind and Msf is the SER. The mass loading 
factor has been fixed to a value = 2 in agreement 
with the CW prescription used in CNll. 


Variable Wind model (VW )— The variable wind model 
(VW) from CNll has introduced a more flexible sub¬ 
grid physics compared to the simple recipe of SH03. In 
this paper, we adopt the 1.5ME wind model described 
in CNll. Briefly, the model assumes that all gas parti¬ 
cles in a given galaxy have the same chance to become 
part of the wind. The probability that the gas particle 
becomes a wind particle is based on the values of /3w, 
galaxy mass and galaxy SER. The main parameters are 
the wind load /3w, defined above, and the wind velocity 
Uw — both have been constrained by observations which 
express these two parameters in terms of the host galaxy 
stellar mass, AA, and the galaxy SER. 

The wind velocity is calculated as a fraction of the 
escape speed from the host galaxy, Uw = C'^esc, where 
C = 1.5 for momentum-driven and C = 1 for energy- 
driven winds in the current setting. The empirical rela¬ 
tion between galaxy SER and Ugsc becomes. 


SER = 1.0 


which is consistent with observations (e.g., [Martin||2QQ5[ 
Weiner et al.||2QQ9 ). The mass loading factor, is as¬ 
sumed to represent the energy-driven wind (/dw oc v~l) 
in the low-density case, n < and the momentum- 

driven wind (/dw cx: u“^) for n > Hence, for the low 

density gas, i.e., away from starforming regions (but in¬ 
side the galaxy!), we apply the parameters of the energy- 
driven wind for the gas particles. This procedure requires 
that simulations compute and Uw using an on-the-fly 
group finder, which, in our runs, is a simplified version 
of the SUBEIND algorithm (jSpringel et al.||2001). 

When an SPH particle is converted to a wind parti¬ 
cle, it receives a Uw kick and decouples from hydrody¬ 
namic forces. The direction of the kicks are chosen to be 
preferentially perpendicular to the angular momentum 
of the particle such that, overall, the winds are launched 
above or below the galactic plane. The wind particles are 
turned back to normal SPH particles when the ambient 
gas density becomes lower than or when they 

travel a distance longer than 20 kpc/h, whichever comes 
first. 
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TABLE 1 

Summary of the different simulation runs and their properties 


Name 

Initial 

conditions 

Wind 

model 

iVeff" 

7? 

J *-zoom 

(h~^ Mpc) 

R- ^ 

-'Unner 

{h~^ Mpc) 

todm'* 

(lOS/l-lMg) 

'UT'gas^ 

(10® /i-IMq) 

(lO^/i-iM©) 

Cgrav^ 

{h~^ kpc) 

■Uw* 

(kms-1) 

CW 

CR 

SH03 

2 X 512^ 

3.5 

2.7 

37.3 

8.88 

4.44 

0.14 

484 

VW 

CR 

CNll 

2 X 5123 

3.5 

2.7 

37.3 

8.88 

4.44 

0.14 

1-1-5'l’esc 

NW 

CR 

- 

2 X 5123 

3.5 

2.7 

37.3 

8.88 

4.44 

0.14 

0 

UCW 

UCR 

SH03 

2 X 5123 

7.0 

4.0 

37.3 

8.88 

4.44 

0.14 

484 

UVW 

UCR 

CNll 

2 X 5123 

7.0 

4.0 

37.3 

8.88 

4.44 

0.14 

1-1.5l^esc 


^effective resolution in the highest refinement (zoom-in) region, ^comoving radius of zoom-in region, ‘^comoving radius of central region 
used to identify halos and galaxies, "^mass resolution in the DM, gas and stellar component, ®comoving gravitational softening length, 

^outfiow velocity used in the galactic wind model 


2.3. Initial conditions 

The initial conditions have been genera ted using; the 
Constrained Realization (CR) method (e.g., Bertschinger 


§ 


1987 Hoffman & Ribak||1991 Romano-Diaz et al. 2i 


and are those used by Roman’d-Diaz et al. (|2011b[ ), be- 
ing downgraded from 2 x 1024^ to 2 x 512T A CR of 
a Gaussian field is a random realization of such a field 
constructed to obey a set of linear constraints imposed 
on the field. The algorithm is exact, involves no itera¬ 
tions and is based on the property that the residual of the 
field from its mean is statistically independent of the ac¬ 
tual mmiericalvahie_of^ (for more details 


see Romano-Diaz et al. (2011a)). The main advantage 
of such method is to bypass the sampling problem of 
highly overdense regions which are rare and thus require 
to simulate a large volume (> Gpc^) of the universe. 

The constraints were imposed on to a grid of 1024^ 
within a cubic box of size 20 h“^Mpc to create a DM halo 
seed of ^ 10^^ h~^ Mq collapsing by z ^ 6, according to 
the top-hat model. We as sume the ACDM cos mology 
with WMAP5 parameters (Dunkley et al. 2009), = 

0.28, CIa = 0.72, fib = 0.045, and fi = 0.701, where h is 
the Hubble constant in units of 100kms“^ Mpc“^. The 
variance ag = 0.817 of the density field convolved with 
the top hat window of radius Sh~^ Mpc“^ was used to 
normalize the power spectrum. The overdensity in the 
GR models corresponds to ^ 5cr, with respect to the 
average density of the universe. 

We also evolved the same parent random realization 
of the density field used to construct the CR runs but 
without any imposed constraints to represent an average 
region of the universe. In order to compare the effect of 
winds in different environments, we have run the uncon¬ 
strained simulation using the CW and VW models. In 
the following, we refer to individual runs by the name 
of the wind model used in the simulation (see Table 
: CW, VW and NW for constrained runs (CR run^, 
and UCW and UVW for their unconstrained counter¬ 
parts (UCR runs). 

The evolution has been followed from z = 199 down 
to z = 6. Our simulations were ran in comoving coordi¬ 
nates and with vacuum boundary conditions. We used 
the multimass approach and have downgraded the nu¬ 
merical resolution outside the central region using 3 dif¬ 
ferent resolution levels in order to speed up the computa¬ 
tional time. The highest refinement region has a radius of 
3.5 h“^Mpc, having an effective resolution of 2 x 512^ in 
DM and SPH particles. The total mass within the com¬ 
putational box is 6.19 X 10^^ h~^ Mq and the mass of the 
high-resolution region is 1.45 x 10^^ h~^ Mq. Within this 
region, we obtain the particle mass of 3.73 x 10® /i-m© 


TABLE 2 

Properties of the halo and galaxy samples 


Name 

(fe-iMo) 

Afh,z=6® 

(li-iMo) 

^g,z=6® 

CW 

10® 

6004 

3 X 10^ 

509 

VW 

lO® 

5768 

3 X 10^ 

425 

NW 

lO® 

5960 

3 X 10^ 

1969 

UCW 

10® 

9135 

3 X 10^ 

245 

UVW 

10® 

8470 

3 X 10^ 

92 


^Minimum halo total mass considered 
^Number of halos at 2 : = 6 
‘^Minimum galaxy stellar mass considered 
•^Number of galaxies at 2 : = 6 

(DM), 8.88 X 10^ ^“^Mq (gas) and 4.44 x 10^ Mq 
(stars). The gravitational softening is Lgrav = 140h“^pc 
(comoving), which is about 20h“^pc in physical units at 
z = 6. 

In order to avoid contamination of heavy particles from 
the outer refinement levels, we have focused on the cen¬ 
tral inner spherical region within radius of 2.7h“^Mpc 
for our CR simulations, and with 4h“^Mpc for the un¬ 
constrained (UCR) simulations. The size of the inner re¬ 
gion in the UCR runs has been chosen by matching the 
cumulative mass profile with the CR models. The prop¬ 
erties of the different simulation runs are summarized in 
Table m 

Finally, the CW model is switched on at the initial 
redshift, while the VW model is run as NW (i.e with¬ 
out outfiows) until z = 12, at which point the winds are 
switched on. This is done in order to speed up the sim¬ 
ulation until the SF starts to rise and the SN feedback 
becomes important. 


2.4. Group finder 


We use the group finding algorithm HOP (Eisenstein 
& Hut |1998|) to identify halos and galaxies. Halos are 
isolated according to a purely particle density criteria — 
the (local) total particle densities, i.e., DM + baryons, 
are calculated with an SPH kernel. A halo is defined 
as the region enclosed within a total iso-density contour 
of AcPcrit(^): whcrc Ac = 80 and Pcrit(^) is the critical 
density at redshift z. Hence, no geometry of the particle 
distribution is assum ed with this approach. Thi s is a 
change compared to |Romano-Dfaz et al.| (|2011b|) who 
used the stand yd pure DlVl definition for the halo, but 
is exactly as in Romano-Diaz et al. (2014). The benefit 
of including the baryonic contribution to the density field 
when identifying halos, while not significantly modifying 
the resulting halo catalog, is to automatically identify the 
baryonic component inside each halo. In order to avoid 
resolution effects, we adopt a minimum mass threshold 
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Fig. 1. — Left: Evolution of the distribution of wind particles as a function of their outflow velocities in the VW run. The vertical dash 
line corresponds to the fixed velocity of 484kms“^ used in the CW model. Right: The median wind velocities as a function of galaxy 
stellar mass in the VW model at 2 ; = 6. The horizontal dashed line shows the velocity used in the CW model (as above). Filled circles 
represent individual galaxies in the VW run, color-coded by their median SFR. The solid line is showing the median relation between 
wind velocity and stellar mass for these galaxies, whereas the dashed and dotted line represents the expected trends for energy-driven and 
momentum-driven winds respectively. 


of Mh = 10^ 

Galaxies arc identified with respect to the baryonic 
density field with the outer boundary corresponding to 
an iso- density contour of Q.Qln,^ 

that of|Roniano-Dfaz et al. (2011b|) but is as inlRomano 


Uiaz et ah ( |2U14p . T'his dehnition ensures inclusion of 


regions which host star-forming gas (n > as well 

as lower density non-star forming gas (n < which 

is roughly b ound to the galaxy. 

Eollowing Romano-Dfaz et al. ( 2Q14[ ), we also exclude 
galaxies below a minimum stellar mass of 3 x 10^ h~^ 
in agreement with recent observations ( Ryan et al.|2Q14 ), 
corresponding to ^ 70 stellar particles, in order to 
verify that our simulations are not affected by numer¬ 
ical resolution, we have compared our GW run with the 


higher-res olution version presented in Romano-Dfaz et 


al. (2014). We found that both runs are in excellent 


agreement in terms of galaxy stellar mass functions, star 
formation rates and metallicities in galaxies, indicating 
that our simulation runs are numerically converged. 


3. RESULTS: EFFECTS OF WINDS AND ENVIRONMENT 
ON HALO AND GALAXY POPULATION 

In this section, we compare the wind models in the 
overdense and average density environments and analyze 
their effects on the mass functions of halos and galaxies 
, as well as on the intrinsic properties of these objects, 
such as gas and stellar fractions, star formation rates and 
on the mass-metallicity relation at 2: > 6. 


3.1. Wind velocities 

Section 2.2 provided the general description of the wind 
models implemented in this work. We start by verifying 
the kinematic properties of these winds, and especially 
of the VW run. In Eigure[^ we plot the resulting veloc¬ 
ities of the wind particles in the VW and GW runs. The 
left panel exhibits the distribution of wind particles as 
a function of their outflow velocity, at different red- 
shifts. The vertical dashed line is the value used in the 
GW model, i.e., 484 km s“^. Even though, overall, the 
total number of particles in the wind is steadily increas¬ 
ing with time, we observe three interesting trends in the 


VW model: the number of wind particles corresponding 
to velocity at the peak of the distribution is increasing 
with redshift until 2: ^ 8 and falls off thereafter; the aver¬ 
age wind velocity increases with time and the tail of the 
distribution extends gradually towards higher velocities. 
Despite this latter increase, we note that, at all redshifts, 
the majority of wind particles in the VW run have lower 
velocity compared to the constant velocity adopted in 
the CW model. 

The right panel of Eigurej^ displays the relation be¬ 
tween the median and stellar mass in galaxies at 2: = 6 
in the VW run. The value of used in the CW model is 
shown as the horizontal dashed line. Individual galaxies 
are represented by filled circles with color scaling based 
on their SER. The solid black line is the median trend 
calculated from these objects. The black dot-dashed and 
dashed lines are the expected trends for the momentum- 
driven and energy-driven winds respectively. They have 
been calculated from eq.[^ assuming a specific SER (i.e., 
sSER) of 2.5 Gyr“^, in accordance with results shown in 
Eigurej^ 

We find the scaling of wind velocities with stellar mass 
to closely follow the energy-driven case for nearly the en¬ 
tire mass range considered, although the scatter increases 
towards low masses. This is not entirely surprising since, 
in the VW model, winds coming from non star-forming 
regions (n < are assumed to be energy-driven. 

Given our definition of a galaxy, which inclu des g as well 


below the star formation threshold (section 2.4|), these 
regions usually dominate the total gas content m a given 
galaxy. Since all gas particles have an equal probability 
to turn into wind, the majority of wind particles emitted 
from a galaxy come from those low-density regions and 
thus produce energy-driven winds, as seen on Eigure[2 
We also find that, except for the most massive systems, 
nearly all galaxies in the VW model have lower median 
wind velocities than the GW case, in agreement with the 
velocity distributions presented in the left panel. The 
most massive objects appear to have stronger winds in 
the VW case than in the CW case. However, since they 
constitute the deepest potential wells, these objects are 
also the ones that are least affected by winds compared 
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Fig. 2. — Snapshots of the mass distribution within the central 1.3/i“^Mpc (comoving) at 2; = 6 in the three different models considered 
here: Constrained Simulations (CR) with Constant Velocity (CW) wind (top row), Variable Velocity (VW) wind (second row), No wind 
(NW) (third row), and Unconstrained simulations with Constant Velocity (UCW) wind (fourth row) and Variable Velocity (UVW) wind 
(bottom row). Each column represents the view for a different component: Dark Matter (left column), gas (middle column) and stars (right 
column). In all panels, the same color scale represents column density in the corresponding component with an arbitrary normalization 
and is used only for visual representation of the mass distribution. 
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to the intermediate and low-mass galaxies. Thus, 
overall, we expect mechanical feedback from winds in 
the VW run to increase with time, peak at 2: ^ 7 — 8 and 
decline afterwards, as shown by the redshift evolution 
in the left panel. At the same time, this feedback is less 
efficient than in the GW run at all times. Given the 
direct coupling between wind and galaxy properties in 
the VW model, this non-monotonic evolution is accom¬ 
panied with a similar behavior in galaxy gas content 
and SER as we show below. In terms of stellar feedback, 
the VW model, therefore, represents an intermediate 
situation compared to the other two cases, NW and GW. 


3.2. Mass distribution and mass functions 

Eigure shows the mass distribution of DM, gas and 
stars in all the simulation runs at 2: = 6 on a scale 
of 1.3h“^Mpc (comoving), corresponding to a physical 
scale of ^ 186h“^kpc. In the same environment (CR 
or UGR), no substantial differences are expected in the 
DM distribution on large scales since they only differ 
due to the physics of galactic outflows implemented. On 
the other hand, we do observe that the distribution of 
baryons differs between the runs on progressively smaller 
spatial scales. 

Eor a more quantitative analysis, we start by ex¬ 
amining the differential HMEs, ^ = dn/dlogMh, i.e., 
the number of halos per unit volume per unit loga¬ 
rithmic total halo mass interval, as shown in Eigure 
(upper frame). As expected, the HMEs of the differ¬ 
ent wind models in a given environment (GR or UGR) 
are essentially identical over the entire mass range of 
lO^h“^M0 — lO^^h“^M0. Slight differences are ob¬ 
served between wind models because our definition of 
a halo i s ba sed on the total mass which includes baryons 
(section |2.4[ ). 

On the other hand, the effect of the overdensity is 
readily visible from this figure. Eirst, the HMEs for 
the three GR models are shifted up with respect to the 
UGR HMEs because of the overdense region sampled by 
these models. Second, the slopes of the three models 
are also somewhat shallower than that of the theoretical 


Sheth-Tormen slope ( 

Sheth & Tormen 

1999), as expected 

around density peakj 
Third, they extend t 
the UGR halos. The 

5 (e.g., Romano-1 

Jiaz et al. 

2011a). 

oward higher masses, compared to 
HMEs in the UGR runs appear to 

be in reasonable agreement with the |Sheth & Tormen 


deviation from the theoretical prediction. 

We have analyzed the possible causes of this excess of 
low-mass halos in the UGR simulations over the Sheth- 
Tormen HME. Eirst, it could follow from the spurious 
detection of gas-dominated blobs which might be identi¬ 
fied as low-mass halos by HOP. We, therefore, have com¬ 
puted the baryon fraction /b = Mbaryons/AfoM for each 
halo in the UGR runs at z=6, but found that /b is always 
below 40%. Second, the difference is still present even 
when we use only the DM component to find and define 
our halos. Third, as a complementary test, we have also 
analyzed the HMEs in a DM-only version of our UCW 


model (run in a larger volume) presented in Romano- 
Dfaz et al.|(|2011a) and found a good agreem ent with the 
theoretical HMFs from Jenkins et al. (2001) over the en¬ 
tire mass range considered here, for this reason, we are 


confident that the departu re between our HM Es and the 
theoretical prediction from |Sheth fc: Tormen| is due to the 
inclusion of baryonic physics, which is expected to have 
the most effect on low-mass halos (in preparation). 

The bottom panel of Eigure shows the galaxy mass 
functions (GMEs) for all models as a function of galaxy 
stellar mass, AU, at different redshifts. Here we observe 
a clear impact of the outflow model on the shape and 
evolution of the GME, reflecting the effect of winds on 
the galaxy growth process. At every redshift, the nor¬ 
malization (amplitude) of the GME in the NW model 
is the highest overall among the GR runs, which means 
that galaxies of a given stellar mass are more abundant 
in the NW model compared to other two CR runs (GW 
and VW). However, since the CR runs have the same 
underlying HMEs, an alternative way to characterize the 
differences between GMEs is to compare their relative 
shift at a given number density 4> (i.e at a given halo 
mass). In this case, we find that the GME in the NW 
model is shifted towards higher stellar masses, meaning 
that, in a given halo, galaxies have produced more stars 
in the NW case because of the absence of feedback from 
winds with respect to the models including outflows (GW 
and VW). 

In comparison, the effect of winds is prominent in the 
evolution of the VW and GW GMEs at different red- 
shifts. At 2: = 11, VW and NW GMEs closely follow each 
other as expected since these two models are identical 
unti l win ds are turned on at 2: ^ 12 in the VW run (sec¬ 
tion [2^). With decreasing redshift, the outflows become 


stronger which cause both VW and GW GMEs to become 
shallower, especially at the low-mass end. At 2: = 6, the 
VW GME is more closely related to the GW one, but 
their slopes differ, which results in a smaller number of 
low-mass (^ 10^h~^MQ) galaxies and a larger popula¬ 
tion of objects at intermediate masses in the VW run. 
The differences mentioned above are no longer visible at 
the massive end (AA ^ 10^^ h~^ M0) of the GMEs where 
all GR runs agree reasonably well. This means that winds 
are inefficient in removing baryons from the deepest po¬ 
tential wells and thus affect mostly the low-mass end of 
the GMEs, in agreement with previous models and obser- 
vations (e .g., Benson et ^ 2003 iDe Young fc Heckman 


19941 |Mac Low V" Eerrara||1999[ iDcannapieco et al.[!^U0I[ 


!||201ip. 

f the G 


Ghoi & Nagamme| 

4'he evolution or the GMEs in the average density re¬ 
gion (UGW and UVW models) appears to resemble their 
wind model counterparts (GW and VW) in the CR runs, 
but shifted towards lower stellar masses. In particular, 
we observe a flattening with time of the low-mass end 
of the UVW GME compared to the UCW case. How¬ 
ever, the main difference between the two environments 
lies in the low-mass end of the GMEs which appear to be 
steeper for UGR than for CR runs. The difference is espe¬ 
cially pronounced when comparing GW and UCW runs 
which have low-mass end slopes of ^ — 2 and ^ —1 re¬ 
spectively. Because UGR and CR runs only differ by the 
presence of the imposed overdensity, it is clear that these 
trends reflect the difference in the underlying HMEs, as 
seen in the top panels of Eigure This is an impor¬ 
tant result concerning galaxy evolution in overdense en¬ 
vironments at high-redshifts which should have implica¬ 
tions for the expected contribution of low-mass galaxies 
to reionization at 2: > 6, since measurements of the faint- 
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log M* [h-^ Mq] 


Fig. 3. — Top: Evolution of the halo mass functions (HMFs) for the three CR wind runs considered here (CW, VW and NW) and for 
the UCR runs (UCW and UVW) at different redshifts. Mh is the total halo mass including DM+baryons. The dashed line corresponds 
to the theoretical Sheth-Tormen HMF. HMFs are insensitive to the outflow model used but show an overproduction of structures in the 
high-density region with respect to the average. Bottom: Evolution of the galaxy mass function (GMF), as a function of their stellar mass, 
for CR and UCR runs. The dashed lines show the trend expected from a power-law relation scaling respectively as M~^ and M~‘^. In both 
top and bottom panels, the mass functions are normalized by the binsize which is set to a constant 0.3 dex in logM. Error bars correspond 
to the expected Poisson Icr deviation in each bin. 


end slope of the UV L F at hi^h-z indicate a steep trend 


with a slope of ^ —2 (jDressler et al.||2015t |Song et al 
2015) in agreement with the trend we nnd in the average 
density regions. 


3.3. Gas fractions: effect of winds 

Next, we analyze the impact of the different galac¬ 
tic wind prescriptions on the gas content of halos and 
galaxies. Figure 0 (top) displays the halo gas fraction, 
/gas,h = ^gas,h/^^ a functiou of rcdshift in four mass 
bins, from Mh = 10^ to 10^^ h~^ Mq with a bin size of 
1 dex. We keep these mass bins constant in time, which 
means that halos are assigned to one of the bin based 
on their total mass at a given redshift. In each panel, 
solid lines represent the median trend of each model, and 
shaded regions enclose the 20 and 80 percentiles of the 
population in each bin. 


For all models, we observe that the evolution of the 
gas content in halos shows a clear dependence with halo 
mass. Low-mass halos start with a gas fraction close to 
the universal baryon fraction, /bar,o = ^b/^m ^ 16 % 
(dashed line), and exhibit a steeper decline with 2: than 
the more massive ones. As halo mass increases, the curve 
flattens, with the highest-mass bin, > 10^^ 
showing little variations in gas fraction, from z = 9 to 
z = 6, although trends are more noisy in this bin because 
of the low number of objects at each 2;. At later times, 
/gas,h is in the range of ^ 5 — 15%, below the universal 
value /bar,o- Note that the highest mass bin does not 
contain any UCR halos since such massive objects have 
not formed yet in the average density region. 

The CR runs (CW, VW and NW) show clear differ¬ 
ences among themselves, reflecting the impact of the out¬ 
flow model on the gas content of halos. On average, at 
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Fig. 4. — Evolution of gas fraction in halos {top) and galaxies 
(bottom) as a function of redshift in different mass bins. Objects 
are assigned to one of the bins based on their current total mass 
(halos) or stellar mass (galaxies) at each redshift. Bin limits are 
indicated in each panel in units of Solid lines are the 

median trends in each bin. Shaded regions show the scatter around 
the median between the 20 and 80 percentiles of the population. 
In the top figure, the dashed horizontal line indicates the universal 
baryon fraction 


each 2:, gas fractions are highest in the NW model for 
lower mass halos, < 10^^ h~^ Mq, and become lower than 
in the GW and VW models for the most massive objects, 
> 10^^ h~^ Mq. Although it seems contradictory at first 
that NW halos have less gas, we find this trend to be ac¬ 
companied by a similar trend in SERs. The SE is more ef¬ 
ficient and starts earlier in the NW model because of the 
absence of a substantial feedback from galactic outflows. 
The SERs are also mass-dependent with very little SE at 
the low-mass end in all models. As a consequence, since 
halos in the NW run are able to better retain their gas 
compared to the GW and VW, more of it gets consumed 
to form stars in higher mass objects which explains the 
above trend. Such a high efficiency in converting the gas 
into stars is one of the reasons that the NW model can 
be ruled out. 

Eigure 0 (bottom) displays the evolution of the gas 
fraction, Jgas,gai = Mgas/(Mgas + AA), in galaxies in four 
mass bins. The lines and shaded areas have the same 
meaning as in the case for halos (top) but now the galax¬ 


ies are binned in terms of their stellar mass at a given 

2:. 

The comparison between the GR models shows a di¬ 
verse population of objects whose gas fraction evolution 
is much more complex than in halos. The NW model 
shows a monotonic decline of /gas,gai with time, except 
in the highest mass bin for 2: > 9. By ^ ^ 6, the gas frac¬ 
tion is decreasing with increasing mass bin. At the same 
time, /gas,gal is independent of mass in the GW model. 
Eor the VW model, the final /gas,gai is also lower for 
higher galaxy masses. The overall evolution of GW and 
VW models, however, differs substantially. Eor almost 
the entire mass range, the gas fraction stays constant for 
2: > 8 — 9, and exhibit a mild decline thereafter. On 
the other hand, /gas,gal in the VW model shows a sharp 
decline for 2: > 9, followed by a sharp increase, and sat¬ 
uration. In the highest mass bin, another sharp decline 
can be clearly observed. 

Overall, by ^ ^ 6, the NW galaxies appear to be most 
gas-poor, followed by VW galaxies and GW ones. Obvi¬ 
ously, this trend correlates with the efficiency of feedback 
in the form of galactic outflows. The weakest feedback 
is that of the NW model. The GW feedback is strong 
but fixed in time, while that of the VW is increasing 
with time, but not monotonically. So averaged over the 
whole galaxy population, the effect of the VW model 
falls roughly between that of the GW and NW models as 
already mentioned in section |3.1[ 


In comparison with Ghoi & Nagamine (2011) galax¬ 


ies, the spread between various models is smaller at 
2 : ~ 5, and /gas,gal -^7gal’ /gas,gal 0.9 at 

log gal r\j 7. Our GR models show smaller gas frac¬ 
tions for all galaxies already at ^ ^ 6, which is the result 
of elevated SER caused by the higher accretion rates and 
gas supply in the overdense region. This trend persists 
despite the overdense region, especially at the high mass 
end . This is also true in comparison with [Thompson et| 


al. (2014) models. The explanation for the gas fraction 


evolution that we observed in the GR models lies in the 
interplay between the SE and mass accreti on histories 
onto galaxies (e.g., Romano-Dfaz et al.|2014 ). As we dis¬ 


cuss later on, the kinetic luminosity of these winds has 
a substantial effect on the temperature of the IGM and 
the halo gas — an effect which lowers accretion rates 
profoundly. 


3.4. Gas fractions: effect of environment 

As evident from Eigurej^ the evolution of gas content 
in GW and UCW halos and galaxies is nearly identical in 
all mass bins, except in the mass bin 10^^ —10^^ h~^ Mq, 
where the most massive objects are found in the UGR 
runs. Since both models use the same outflow prescrip¬ 
tion, this strongly suggests that, for a given wind model, 
the overdensity has little effect on the amount of gas in¬ 
side these halos and galaxies. 

When comparing VW and UVW halos, we find that 
their gas fractions tend to follow similar trends in each 
mass bin, but with larger relative differences than be¬ 
tween GW and UGW runs. This indicates that the VW 
model couples galaxy evolution with the large-scale en¬ 
vironment more strongly than the GW case. 

However, addressing the effect of the environment on 
gas fraction is difficult from a direct comparison between 
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Fig. 5.— Top: Evolution of gas fraction in galaxies in the CR and 
UCR runs as a function of redshift, similar to Fig. but objects 
are now binned by density contrast S = p/piy — 1 vmere pb is the 
universal mean baryonic density. The bins are defined as: (5 < 0 
for underdense (left panel), (5 > 0 for overdense (middle panel) 
and all galaxies (right panel). Bottom: Mass-weighted average gas 
fractions in galaxies as function of stellar mass at z = 10 (left 
panel), z = 8 (middle panel) and z = 6 (left panel). 

CR and UCR evolution from Fig. [^because objects in a 
same mass bin can reside in different environments. For 
example, low-mass galaxies can be found in underdense 
regions in CR and UCR models, but also as satellites 
around the massive halos in the overdense region. This 
effectively causes mixing among objects of similar masses 
and can thus wash out the differences induced by the 
presence of the overdensity. 

To investigate further the effects of the environment 
on the gas content in galaxies, it is thus easier to di¬ 
rectly divide objects depending on their large-scale envi¬ 
ronment. Figure^ (top) shows the redshift evolution of 
gas fractions in gmaxies in different density bins for the 
CR and UCR models that include outflows. The bins are 
defined based on the smoothed relative density contrast 
S = p/pb — 1 where pb = ^bpcrit is the universal mean 
baryonic density. The smoothed galaxy density field p 
has been obtained by convolving the local galaxy densi¬ 
ties calculated using an SPH method with a top-hat filter 
of 250 h~^ kpc scale. Note that we use the universal mean 
baryonic density pb in our definition of 5 instead of the 
mean baryonic density inside the simulation box in or¬ 
der to have a common normalization factor for both CR 
and UCR models. Galaxies are divided between under¬ 
dense {6 < 0, top left panel) and overdense ((5 > 0, top 
middle panel) bins and we also show the evolution for 
the entire galaxy population for comparison (top right 
panel). Based on our galaxy density definition, we found 
no galaxies with ^ > 0 in the UCR models. The compari¬ 
son between CR and UCR runs is therefore only possible 
for galaxies residing in underdense regions. 

The top left panel of Figure shows that the evolution 


of gas fraction in these galaxies residing in underdense 
regions follows a close trend in both CW and UCW mod¬ 
els with a modest decline from ^ 65-70% at z = 12 to 
^ 55% at 2 ; = 6. Differences between the two models are 
not significant given the noise introduced by the small 
number of galaxies formed in the UCW model. In com¬ 
parison, the gas fraction in galaxies residing in overdense 
regions in the CW model exhibits a much steeper decline 
from z = 10 to z = 6 (top middle panel) due to the ac¬ 
celerated evolution in these regions. The UVW galaxies 
show a more complicated behavior, but decline as well 
after 2 : ^ 8, in tandem with CW. The evolution for the 
entire galaxy population (top right panel) shows larger 
differences between CW and UCW due to the mixing 
of both underdense and overdense environments in the 
CW case but are still showing a similar decreasing trend 
with redshift and comparable gas fraction levels at all 
The VW and UVW models in the underdense bin ex¬ 
hibit much less correlation than CW and UCW among 
themselves, as already noted in the previous section. 

The similar gas fraction evolution for CW and UCW 
galaxies residing in the same environment further con¬ 
firms that the gas content is mainly affected by the dif¬ 
ferent wind models and not by the presence of the over¬ 
density. The gas fraction in VW and UVW models dif¬ 
fer more among themselves, which can be explained by 
weaker winds they produce. To more fully understand 
the reason for this effect, we need to investigate how the 
galaxy gas fractions depend on the outflow prescription 
at a given redshift. Figure (bottom) shows the mass- 
weighted average galaxy gas fractions as a function of 
galaxy stellar mass in all the simulation runs at ^ = 10, 
8, and 6. Since we are only interested in the average 
trends, the curves have been smoothed with a boxcar fil¬ 
ter of size 0.5 dex in stellar mass to reduce Poisson noise 
caused by the small number of objects in each mass bin. 

We find significant differences among the various wind 
models at all ^ but the main discrepancy can be ob¬ 
served between CW and the other two runs, VW and 
NW. In particular, there is a clear trend of decreasing 
gas fraction with increasing stellar mass in both VW 
and NW, which is not seen in CW. The latter shows 
nearly constant gas fractions across the entire mass range 
= lO^’^-lO^^ h~^ Mq. The UCW model also exhibits 
nearly mass-independent gas fractions at all z, which fol¬ 
lows the CW trend as expected. This is in agreement 
with results shown in the top panels of Figure and 
Figure although there are clear signs of divergence 
between the two models for ^ < 8. These results thus 
suggest that the likely explanation for the apparent sim¬ 
ilarities in galaxy gas content between CW and UCW 
can be attributed, for the most part, to the constant 
and strong outflow prescription used in both runs. Since 
winds scale independently of galaxy properties (mass and 
SFR) in this outflow model, the gas content inside galax¬ 
ies varies weakly with stellar mass. As a consequence, the 
CW and UCW runs produce galaxies with nearly con¬ 
stant gas fractions, which explains the similar behavior 
observed between the two models despite their different 
environments. 

3.5. Star formation histories 

To construct the star formation history (SFH) for a 
given galaxy selected at 2 ; = 6, we identify the pro- 
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Fig. 6. — Star formation histories for galaxies in different stellar mass bins. Galaxies are selected based on their stellar mass at 

z = 6 (which is different than selecting galaxies based on their mass at every redshift) and their ancestors are identified at earlier redshifts 
to be the most massive object among all the progenitors. The bin limits are indicated in each panel in units of h~^ Mq. Solid lines show 
the median SFR evolution of all the ancestors at each 2 ;. The shaded regions are bounded by the 20 and 80 percentiles. 


gcnitors at earlier redshifts by matching the particle 
unique identification number (IDs). The ancestor galaxy 
is taken to be the progenitor containing the largest num¬ 
ber of particles. We repeat this process iteratively to 
obtain the list of ancestors of the selected galaxy, from 
2 ; = 6 to ^ = 12. Eigure shows the redshift evolu¬ 
tion of galaxy SERs in different mass bins. Galaxies are 
assigned to one of the bins based on their stellar mass 
at z = 6, and their SEH is then followed back to 
2 : = 12 using the above method. We stress that this mass 
binning is different from the one in which galaxies are se¬ 
lect ed b y mass at each redshift, as used in sections |3.3| 
and |3.6[ Eor this reason, attempting to correlate SER 
and gas fraction evolution cannot be performed using a 
straightforward comparison between figures and as 
they represent different galaxy populations at each red¬ 
shift. 

We find a strong dependence of SEHs on stellar mass in 
all models. Overall, massive galaxies exhibit a strongly 
rising SER, while the lowest mass galaxies maintain fiat 
SERs on average. These trends are well correlated with 
the mass accretion rates onto the galaxies, analyzed in 
Romano-Dfaz et al. (2014). More specifically, in the low¬ 
est mass range, 7.5 < log < 8.5, SERs are es¬ 

sentially constant with z for all models, and lie in the 
range of ^ 0.1 — 1 Mq yr“^. The intermediate mass bins, 
8.5 < log < 9.5 and 9.5 < log < 10.5, ex¬ 

hibits an increase in SER by a factor of ^ 7 — 8 between 
2: = 12 and 2: = 6. Einally, the most massive galaxies dis¬ 
play the strongest increase, roughly a factor of 30 over 
this redshift interval. Note, however, that small SE levels 
seen for low-mass galaxies likely result from poorly re¬ 


solved outflows, with only few particles for these objects. 
Eor this reason, we mostly focus on SERs in intermediate 
and high-mass galaxies in the following. 

At high 2 :, the SERs for the most massive galaxies ap¬ 
pear similar among all the GR runs (Note, these galaxies 
are absent in the UGR runs). The rates start to diverge 
below 2 : ^ 9, with the GW model becoming the domi¬ 
nant one, followed by the VW and the NW, each lower 
by a factor of 2 - 3 by 2 : = 6. Eor these galaxies, the 
SER has reached ^ 10^ Mq yr“^, bringing the m into the 
ob served regime of LBGs at high- 2 : (see also Yajima et 


al. 2015). We also observe a similar trend in the inter 
mediate mass bins, where SER in the NW model start 
to decrease after 2 : ^ 8, while models including galactic 
outflows continue to increase their SERs. By the end of 
the simulations, the GW model has the highest median 
SER with ^ 30 — 40 Mq yr“^ for 9.5 < log Mf < 10.5, 
and - 600 - SOOMQyr-^ for 10.5 < log < 11.5. 

The corresponding median SERs in the NW run are 
^ 10MQyr“^ and ^ lOOMQyr”^, and the VW model 
ends up with intermediate values between the GW and 
NW runs. 

These results indicate that, even at high-z, galactic 
outflows do play an important role in controlling galaxy 
growth in over dense regions. In the NW run, the ab¬ 
sence of these outflows means the absence of a strong 
feedback which can delay or quench star formation. As a 
consequence, the peak of star formation in the NW run 
happens at higher redshifts, 2 : > 6, as noted above, and 
thus the available gas supply for SE in galaxies is already 
consumed by 2 : = 6. Subsequent gas accretion could in 
principle increase SERs again at later times in the NW 
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Fig. 7. — Me dian specific star fo rmation r ate (sSFR) as a func tion of stella r mass deft) and reds hift (right). Observational data points 
are taken from|Labbe et al.|(|2013|) at 2 : = 8,|Salmon et al.|(|2015| at 2 ; = 6 and ' ^ 

frame) is an observational constraint at 2 ; ~ 5 — 5.5 trom the V iJDS survey ^ 


[ Duncan et al. 

(|2014|) at 2 ; = 6 — 7. The yellow hexagon (left 

lasca et al.j (jz 

LUbj), see their Fig. 3). The sSFR-stellar mass 


[Salmon et al.| ( [2Q15| >. In the right panel, the solid lines correspond to the median sSFR of all the galaxies at each redshift. The dashed 
lines are the median sSFR for galaxies with stellar masses in the range 10^ — 10^'^ h~^ Mq which approximately matches the mass range 
of observed galaxies corresponding to the data points. 


run but, at ^ = 6, we see that accretion has not yet been 
able to replenish the gas content in galaxies to maintain 
high SFRs. For this reason, SFRs are higher in CW and 
VW models at z = 6 since the feedback from winds in 
those models acts to effectively delay the peak of star 
formation activity. 


This is in agreemen^t with|Choi fc Nagamine (2011) and 


Angles-Alcazar et al.| (|2014|) which have analyzed galaxy 
evolution in average density regions, but lack massive 
galaxies at these redshifts. However, since galaxy evo¬ 
lution is accelerated in the overdense regions, the peak 
of SF activity happens at earlier times than in average 
density environment. When comparing CW and VW 
models, results published in the literature show that the 
CW SFRs are generally the lowest o nes, as measured in 
simul ations at low 2 ; ^ 3 — 4 (e.g., |Choi & Nagamine 
201 Ij). However, the cosm ic SFR evolution presented in 
Choi V Nagamine (2011) (see their Fig. 5) shows that 
the ditterence between CW and VW SFRs decreases with 
increasing redshift at 2 : > 6, and, that these two models 
eventually reach similar SF levels at high- 2 ;. 

3.6. Specific star formation rates 

For a more detailed picture of the SF efficiency in dif¬ 
ferent runs, we investigate the specific star formation rate 
(sSFR), which represents SFR per unit mass, and essen¬ 
tially indicates the production rate of stellar mass per 
stellar mass. Figure displays the median sSFR-stellar 
mass relation in galaxies at 2 ; = 6 (left frame) and the 
sSFR evolution with red shift (right frame) in all runs. 
Observational data from Labbe et al. ] (imi 
Duncan et ^ (|2014|) (z = 6 a nd 7),ISalnion et al.| (|2015|) 
\z = 6) and Tasca et al.|(2015|) {z = o.5) have been added 
to both panels. In the left panel, galaxies are binned by 
stellar mass with 0.3 dex interval. 

We observe a nearly constant sSFR with in the NW 
and VW models, with sSFR ^1.5 and 2.5 Gyr“^, respec¬ 
tively. At the same time, the CW model shows an overall 
positive correlation between sSFR and M^. This upward 
turn above the characteristic mass of r\j 10 ^ h ^ Mq is 
in contradiction with observational trend shown on Fig¬ 
ure [71 


We also find that the sSFR in the CW run fall within 
the observed range, but the exsistense of the peak in the 
mass range of ^ 10 ^ — 10^^ h~^ does not agree with 
the observations of [Salmon et al.j (2015|). On the other 
hand, for the NW run, the overall trend roughly agrees 
with the observed one, but the sSFR lies outside of the 
observed range. The best fit to observations is that of 
the VW model — both the trend with and the sSFR 
range agree, although the median observed fit is slightly 
shifted away from the model values. 

The right-hand panel of Figureshows the evolution of 
the sSFR with redshift. Solid lines correspond to the me¬ 
dian trends for the whole galaxy population. The dashed 
lines show the evolution only for galaxies with stellar 
masses in the range 10^ — corresponding 

approximately to the estimated masses of objects tar¬ 
geted by observers. We find the median evolution of all 
galaxies to be dominated by the low and intermediate- 
mass objects because of a higher number of these galaxies 
compared to the massive ones. The median sSFR of all 
galaxies in the CW and UCW runs are in good agreement 
— again underlying that gas content and SF in low-mass 
galaxies are not affected significantly by the overdensity. 

The sSFR in the NW run is the highest for z > 7, 
where it falls below the VW model. In all models except 
VW, the sSFR drops by almost an order of magnitude, 
from ^ 10Gyr“^ at 2 : = 12 to ^ 1 Gyr“^ at ^ = 6. The 
sSFR evolution in the VW case exhibits a trend similar 
to the one found for the gas fraction in low-mass galaxies 
(Fig.0. As we have argued in the previous section. This 
trencTis caused by scaling of the feedback strength with 
the SFR in the VW run. At high- 2 :, outflows quench the 
SF, which in turn decreases the feedback from the winds 
and allows an elevated SFR at later times. The median 
sSFR evolution for the whole galaxy samples does not 
reproduce correctly the observations for 2 ; < 8. The ex¬ 
planation of the mismatch lies in the dominant contribu¬ 
tion from low-mass galaxies (M^ < IO^/i^^Mq) to the 
global sSFR. When comparing the evolution for galax¬ 
ies in the observed stellar mass bin (dashed lines), both 
GW and VW, as well as UGR, runs appear to be in a 
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Fig. 8.— Left: Gas (top) and stellar (bottom) metallicities of halos (including gas in galaxies) as a function of halo total mass at 2 ; = 6 
in all simulation runs. Right: Gas (top) and stellar (bottom) metallicities of galaxies as a function of stellar mass at the same redshift. The 
solid lines indicate the median values, w hile the colored regions correspond to the values between the 20 and 80 percentiles. Metallicity 
profiles of the IGM gas are shown in Fig. 


good agreement with the observed increase of sSER with 
redshift. 


3.7. Mass-metallicity relation 

Next, we focus on the mass-metallicity relation in halos 
and galaxies. Eigure (left frames) shows the metallic- 
ity of gas and stars m halos as a function of the total 
halo mass at z = 6. The highest metallicities are found 
in the most massive halos because their star formation 
starts earlier and they exhibit a higher SER compared to 
lower mass objects. We observe a rapid increase in the 
metallicity in the low-mass halos, and its leveling off at 
the intermediate/high mass end. The GW and VW runs 
reach close to solar metallicity for Mh ^ 10^^ h~^ Mq 
by the end of the simulations. Eor smaller Mh, the GW 
and VW models display a slow decline in Z, which steep¬ 
ens dramatically below Mh r\j h-^ M© for VW, and 
below ~ 10^ h~^ M© for the GW. This ‘knee’ is moving 
gradually toward higher masses with redshift. 

While halos with Mh ^ 10^^ — 10^^ h~^ M© have a me¬ 
dian ZjZ^ ^ 0.6 (gas) and ^ 0.8 (stars), the lower mass 
halos, with ^ 10^ — 10^ h~^ M©, show ZjZ^ ^ 0.01 — 0.03 
(gas) and ^ 0.03 (stars), respectively. Note that the NW 
model achieves the highest metallicity in all Mh, for both 
gas and stars. Eurthermore, compared with the halo gas 
metallicities, stellar metallicities display higher median Z 
— this means that the gas metallicity is dilut ed by the 
cold accretion fro m the IGM, in agreement with |Romano- 
Dfaz et al. (2014|). 

Galaxy metallicities in the gaseous and stellar compo¬ 
nents are shown in Eigurej^ (right frames). As expected, 
the NW model metallicity exceeds substantially those of 
GW and VW in both components, due to the metals be¬ 


ing locked in galaxies and their immediate vicinity. A 
possible exception may be the highest mass bin for the 
ISM, where GW model shows marginally higher metallic- 
ity. The bimodal behavior observed for the halos, can be 
seen here as well. Stars also remain more metal-rich than 
the gas, again in tandem with the halos. Euture obser¬ 
vations will constraint these models. We point out that, 
in Eigure we have only computed the average metal- 
licity inside individual galaxies and ignored the scatter 
in metallicity in a given galaxy. 

4. RESULTS: EFFEGT OF WINDS ON ENVIRONMENT 

In this section we focus on the effects of winds on the 
thermal evolution and metal enrichment in the circum- 
galactic and intergalactic mediu m. W ind effects on the 
ISM has been discussed in section l3^ 

4.1. Temperature and metallieity maps 

The temperature maps of the IGM gas, i.e., gas out¬ 
side galaxies and their host halos, in the inner central 
region of the simulation box clearly display the growing 
difference between the wind models with redshift. Eig¬ 
ure (left frames) underlines profound differences exist¬ 
ing oetween the GR runs already at ^ ^ 10. The NW 
run lacks gas with T > few x 10^ K, except for a few 
spots near the central density peak corresponding to the 
growing massive halo. The warmest gas can be found in 
the large-scale filaments, and it roughly delineates them 
in all GR models. The few x 10^ K gas is more closely 
associated with large-scale structures in the VW case, 
because the feedback from outflows is lower compared to 
the GW run, where this gas “spills out” of the filaments. 
The rest of the volume hosts gas with T < 10^ K. 
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Fig. 9.— Temperature (left) and metallicity (right) maps of the IGM gas (excluding the gas inside halos) within the central 2.7 h ^Mpc^ 
inner region, at 2 ; = 10, 8 and 6 (top to bottom) for the three CR models (CW, VW and NW). 


ucw uvw 



Fig. 10. — Temperature (left) and metallicity (right) maps of 
the IGM in the UCR runs at 2 ; = 6 within the central 4.0 h“^Mpc^ 
inner region. The color scales are the same as in Fig. 

The redshift evolution of these maps shows gradually 
heating and cooling regions, of few x 10^ K gas at z = 8, 
and few x 10^ K gas at z = 6. By the end of the run, 
the massive central halo is surrounded by a hot bubble 
of ^ 10^ K gas in CW and VW wind models. Away 
from the most massive halo, the filament temperature 
drops, most profoundly in the NW model. By ^ ^ 6, 
the CW gas in nearly the whole volume is heated up to 
10^ — 10^ K. In this model, the 10^ — 10^ K gas is present 
in the main filament, and ~ 10^ K gas in other filaments. 
In VW and NW models, the gas outside the filaments 
remains cold. 

Hence, at higher redshifts, z ^ 10, we observe increas¬ 
ing temperatures of the IGM gas, along the sequence 
from NW, to VW, and to CW. With decreasing z, the 


central bubble heats up, and so do the filaments. The 
outflows in the CW model are capable of heating up the 
IGM, while the VW outflows have a dramatically lower 
impact on the IGM temperature at these redshifts, but 
nevertheless heat it up toward z ^ 6. 

Figure (right frames) shows maps of the metallicity 
distribution in the IGM gas at similar redshifts, for all 
the CR models. At ^ = 6, the IGM peak metallicity 
has reached ^ 0.1 Zq inside the central hot bubble and 
along the main filament corresponding to the regions sur¬ 
rounding the massive central structure in the CW and 
VW models. The peak metallicity in the NW is a fac¬ 
tor of a few lower even though the SFRs are higher in 
this model at high-^. The extended region in the CW 
model shows a negative metallicity gradient away from 
the central bubble and the main filament. Still, most of 
the volume is polluted by metals and has Z > 10' -2^0. 
The VW exhibits a much lower metallicity outside the 
filaments, Z < 10' whereas the NW IGM remains 
mostly pristine. The reason for this is the absence of 
outflows, which causes metals to be locked in galaxies 
and their halos. In comparison, solar metallicities have 
been reached only within galaxies of the NW model (see 
Fig.|^. 

Figure displays the temperature and metallicity 
maps of the IGM in the UCR runs at z = 6, and can 
be compared directly with the CW and VW maps in 
Figure!^ as both figures use the same color scales. Since 
the UCR runs represent an average density region and 
lacks massive halos (e.g.. Fig. |^, both the peak and av¬ 
erage IGM temperatures are dramatically lower than in 
the CR runs. The peak metallicity in the IGM and its 
average are also lower compared to the CR runs. Thus 
the impact of the overdensity is readily visible from these 
two figures — the average density universe remains rela¬ 
tively cool at z ^ 6. 
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Fig. 11.— Temperature (top panel) and metallicity (bottom 
panel) radial profiles of the IGM gas at 2 : = 6 calculated from the 
center of the computational domain to the radius i^inner of the 
inner high-resolution region (see table Q. 

4.2. Temperature and metallieity profiles 

A more quantitative comparison between the differ¬ 
ent runs on the state of the IGM at ^ = 6 is dis¬ 
played in Eignre pT] which shows the spherically-averaged 
IGM temperature (top) and metallicity (bottom) pro¬ 
files from the center of the computational box to the 
radius Rinner of the inner high-resolution region, where 
dinner = 2.7 h“^Mpc and Rinner = 4h“^Mpc for the GR 
and UGR models respectively. The profiles are computed 
by averaging the temperature and metallicity of the IGM 
gas on spherical shells of thickness 50h“^kpc. 

Both the temperature and metallicity profiles agree 
with the qualitative results from Eignres [9| and \To\ The 
GR runs have similar IGM temperature within the cen¬ 
tral lh“^Mpc, corresponding approximately to the size 
of the imposed constraint which is visible as a hot bubble 
on Eigurej^ In this central region, the IGM temperature 
reaches T ^ 10^ K which is consistent with the virial tem¬ 
perature for a halo of mass ^ 10^^ Mq, corresponding to 
the mass of the central halo that forms in the overdense 
region. On the other hand, the IGM temperature in the 
UGR runs remains constant with radius at much lower 
values, around ^ 10^-^ - 10^ K and ^ 10^-^ - 10^ K for the 
UVW and UCW models respectively. Thus, it appears 


that shock-heating caused by the presence of the massive 
structure in the overdense region is largely responsible 
for the elevated IGM temperature within 1 h“^Mpc in 
the GR runs. This also explains why the GR runs show 
similar IGM temperature in the central region as the dif¬ 
ferences in temperature due to the different wind models 
(which are much smaller) get washed out. 

The metallicity profiles (Eigurej^ bottom) are in gen¬ 
eral much flatter than the temperature ones. The lowest 
metallicities are seen in the NW run, followed by the VW 
and GW models, in agreement with Eigurej^ The UGW 
IGM metallicity is smaller than both GR models includ¬ 
ing winds (GW and VW) within ^ 2h“^Mpc and then 
rises on larger scales to reach levels comparable with the 
GW run. This increase is associated with the presence 
of a massive filament in the UGR runs (Eigure 10), as a 
similar rise can be seen in the UGW temperature profile 
at the same distance from the center. 


4.3. Gas phases at z = Q 

A supplementary global view on the temperature and 
metallicity of the gas can be taken from Eigure where 
we show the gas phase diagrams in the temperature- 
density plane at 2 : = 6 for all gas particles found in 
the high-resolution region. Gas particles are also colored 
based on their metallicity. The densities are now normal¬ 
ized by the average baryon density in the high-resolution 
region, as opp osed to using the universal value like we 
did in section 13.41 This normalization is thus different 
between GR and UGR runs since the average density is 
higher in the overdense region. Eor this reason, the verti¬ 
cal dashed line, which separates the non-starforming (to 
its left) from starforming (to its right) gas, is slightly off¬ 
set in the bottom panels (UGR) compared to the one in 
the top panels (GR). We observe the highest metallici¬ 
ties, around solar values, in the starforming gas which 
is confined deep within galaxies. One can easily identify 
the wind gas in the lower right part of the diagrams, as 
it is present in the VW, GW, UGW and UVW cases, and 
is metal-rich, but not in the NW case. 

The fraction of low-density gas polluted with met¬ 
als, (green and red) compared to pristine (blue) and 
the highest metallicities attained there, depend strongly 
on the wind model. The metal-rich, low-density gas 
originates in the galactic outflows, and heats-up by the 
shock-stopping of these high-speed outflows. This gas 
lies at temperatures of ^ 10^ — 10^ K gas and densities 
logPb/Pb ^ 4. Here the baryon density pb is normalized 
by the average density in the high-resolution box. We 
note that a large amount of the ^ 10^ K starforming gas 
found in our models at high densities is due to the de¬ 
layed star formation algorithm, the Pressure model, used 
(section 2.1), and because of the absence of H 2 cooling in 
our simulations. The amount of the heated gas increases 
substantially from the NW, to VW and GW runs. The 
hot and low density pristine gas at T > 3 x 10^ K in the 
GR runs consists mostly of virially shock-heated mate¬ 
rial infalling onto the most massive halo, and, as such, is 
not observed in the UGR models. The amount of inter¬ 
mediate metallicity, ZjZ^ ^ 10“^, gas at T > 10^ K, is 
steadily increasing in all GR models. 

Eigure [^provides the metallicity distribution for dif¬ 
ferent wind models in various gas phases — the IGM, 
halo gas, and the ISM, as well as the total distribu- 
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Fig. 12.— Distributions of all gas particles found in the high-resolution region in the temperature-density plane at 2 : = 6 in all the 
simulation runs: (a) CW, (h) VW, (c) NW, (d) UCW, and (e) UVW. Gas particles are colored based on their metallicity as indicated by 
the color scale on the right. The density is normalized by the average baryon density in the high-resolution computational box which is 
different between CR and UCR runs. For this reason, the vertical dashed line corresponding to the star formation threshold density 
is offset in the bottom panels (UCR runs) compared to the one in the top panels (CR runs). 
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Fig. 13. — Metallicity distributions for different gas phases at 
z = 6 in the high-resolution region: All gas particles (top left), 
IGM (top right), halo gas (bottom left) and ISM (bottom right) 
for all the simulation runs: CW (red), VW (green), NW (blue), 
UCW (black) and UVW (grey). 


tion for gas particles in the high-resolution region. The 
ISM includes all the gas inside galaxies, the halo gas 
includes the gas within halos but outside galaxies, and 
the IGM includes the gas outside the halos. In the 
IGM panel, we observe that the GW models for CR 
and UCR exhibit the highest metallicities, basically in¬ 
distinguishable from each other, with maxima located at 
\og{Z/ZQ) ^ (—0.8) — (—0.7). The least efficient injec¬ 


tion of metals happens in the NW model (the maximum 
at \og{Z/ZQ) ^ —2.2), while VW occupies broadly the 
region between the CW and NW models. It displays a 
nearly flat top at \og{Z/ZQ) ^ (—2) — ( — 1). On the 
other extreme, the metal distributions in the ISM are 
basically identical, and the maxima are found between 
\og{Z! Zq) ^ ( — 1) — 0. The metal distributions in the 
halo gas are flat-top for NW and VW models, and show 
maxima at \og{Z/ZQ) ^ —1 for the UCW, and —0.7 
for CW. From Figure we conclude that the ISM is 
polluted by metals in all models. Moreover, there is a 
large difference in expelling the metals to the halo and 
the IGM. The CW and UCW models appear to be the 
most efficient in this process in a given environment. 


4.4. Gas phases evolution 

To take an alternative view at the gas evolution, we 
subdivide the gas phases in a different way. In this 
section, we alter definitions for the IGM and other gas 
phases used in the previous section, in order to facilitate 
the comparison with previous works. Figure com¬ 
pares the evolution of various gas phases, namel 3 ^iffuse, 
warm-hot-hot, IGM and condensed gas. These phases 
are defined based on the temperature of the gas and the 
density fluctuation S = ph/ph — 1. For th i s pur pose, 
we employ the definition s from [Dave et al.J ( |1999| ) and 
Choi & Nagamine (2 Q11[ ). The ‘dittuse' gas is defined by 
T < 10^ K and 6 < 1000, the ‘hot’ gas by T > 10^ K, and 
the ‘warm-hot’ gas by 10^ <T < 10^ K. The ‘condensed’ 
phase consists of non-star forming gas with T < 10^ K 
and 6 > 1000, and of a starforming gas as well as stars. 

We observe that the mass fraction of the ‘condensed’ 
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Fig. 14.— Evolution of the mass fractions of the four different baryon phases (hot, warm-hot, diffuse, and condensed — see section [4^ 
for definitions) for the CR and UCR simulations (see the text for definitions). Panel (b) shows the addition of the two phases, ‘warm-hot + 
hot’, and panel (c) shows the total mass in the IGM, i.e., ‘diffuse -h warm-hot -h hot’. Most of the UCR baryons are in the ‘diffuse’ phase, 
while in the CW the majority is found in the ‘warm-hot’ phase. The ‘condensed’ phase is also more populated in the CR runs, especially 
for the NW model. 


phase is increasing with time for all models, especially 
for the NW, while the IGM (non-starforming) fraction 
decreases. The fraction of a ‘diffuse’ filamentary gas is 
decreasing faster for the GW model than for the VW one. 
At the same time, both the UGW and UVW diffuse gas 
fractions remain very high at all 2 :, while the condensed 
fraction stays very low. While in the warm-hot-hot phase 
the NW and VW evolution are similar, in the condensed 
and IGM phases it is the VW and GW runs which evolve 

in tandem. _ 

While the ma ss fractions in our UGRs and IGhoi fcl 


Nagamine (2011) models appear to be reasonably similar, 
the CK models differ substantially, both from the UCR 
cases and among themselves. This difference between 
the CR and the UCR models is expected, because of the 
overdensity present in the CR runs. But we also observe 
that the NW model exhibits a much higher fraction of 
cold gas, as it evolves toward 2 : ^ 6, and consequently, 
exhibits smaller fraction of IGM gas. The mass fraction 
of the condensed phase, which also includes stars, in the 
NW run reaches ^ 50%, while for the VW and CW mod¬ 
els it is only ^ 32%. The diffuse component shows that 
the fraction drops to ^ 43% in NW, to ^ 51% in CW 
and to only ^ 57% in VW. 


5. DISCUSSION AND CONCLUSIONS 

We have used high-resolution cosmological zoom-in 
simulations of average and overdense regions in the uni¬ 
verse in order to follow their halo and galaxy population 
evolution at z >6. Specifically, we compared this evolu¬ 
tion applying different galactic wind models in order to 
analyze the effects of outflows on: 

• galaxy evolution, including their environment. 

• thermodynamic state of the IGM. 

• the spread of metals in the ISM, DM halos and the 

IGM. 


We have tested three different prescriptions for galac- 
tic winds: a const ant-velocit y wind (CW) model from 


Springel & Hernquist 

(2003 

, a variable-velocity wind 

(VW) model based on the 

Choi & Nagamine] ( 

2011) 

method and a model without outflows (JN W). Each model 


has been followed down to 2 : = 6. In order to get an 
insight into environmental dependence of galaxy evolu¬ 
tion, the CW prescription for galactic outflows has been 
applied to both a constrained (CR) and unconstrained 
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(UCR) simulations, representing overdense and average 
regions in the universe. 

The most general conclusions which emerge from this 
work are that, regardless of the applied wind prescrip¬ 
tion, (1) galaxies in the overdense region evolve much 
faster than their av erage-region counterparts (e.g., Fig¬ 
ures]^ & see also |Romano-Dfaz et al.]|2014|), and (2) 
the Kw-mass end of the galaxy mass function is shal¬ 
lower in the overdense regions — an important issue to 
be addressed elsewhere. This is a direct corollary of the 
imposed constraint which has been designed to collapse 
by 6, based on the top-hat model. As a consequence, 
we observe an accelerated evolution of its environment 
withi n the correlation length of the const rained halo 
(e.g.,|van de Weygaert & Bertschinger||1996|). DM halos 
in the vicinity of the constraint experience a faster evolu¬ 
tion compared to their average-density counterparts, in¬ 
ducing earlier baryon collapse, setting the stage for star 
formation. We find that the galaxy population in such 
regions is much more evolved by the end of simulation. 
This is indicated not only by the presence of more mas¬ 
sive galaxies (up to 1 — 2 x 10^^ h~^ Mq in stellar mass) 
which are absent in the average density region (where 
but also by a substantially larger 
number of galaxies over the mass range considered here. 

Finally, we observe that galactic winds have a substan¬ 
tial effect on galactic environment, which includes mod¬ 
ifying its thermodynamic state and metallicity. While 
CWs show little dependence on the overdensity, the 
corollaries of VWs are more diverse. This is clearly re¬ 
lated to the amount of energy deposited by the VWs in 
DM halos and IGM, which varies with galaxy proper¬ 
ties, and, in general, is smaller than that deposited by 
the CWs. In other words, the VWs appear less “destruc¬ 
tive” than their constant velocity counterparts, and more 
fine-tuned to the galaxy growth mode in the form of cold 
accretion flows — their overall effect differs between over- 
dense and average regions. On the other hand, the CWs 
weaken the cosmological filament flows around galaxies, 
irrespective of surrounding densities. 

Galactic winds, as a form of feedback, serve as main 
agents in distributing recycled, metal-enriched gas, which 
is used to form new stars. The effect of winds onto the 
evolution of our modeled galaxies can be most clearly 
observed in the galaxy mass function (Fig.[^ bottom 
panelT gas fractions (Fig-i and star formation rates 
(Fig.^. The IGM properties can be affected by galactic 
winds, so comparing the observed IGM properties would 
provide the clearest justification for the wind models. 
However, the IGM at z > 6 is not fully reionized, and the 
conventional IGM observation technique, e.g. quasar ab¬ 
sorption line measurement, cannot be used in this epoch. 
This limits the IGM observation resources that can con¬ 
straint our wind models. However, the pre-reionized HI 
status at z > 6 provides an interesting implication that 
can discriminate between our GW, VW and NW models, 
as discussed below. 

The model with no-wind feedback produces the largest 
number of galaxies because the gas is hardly expelled 
from their DM halos. In this case, the feedback is lim¬ 
ited to the thermal feedback by the SN, and, therefore, 
contributes to the gas and galaxy survivals within the 
DM substructure. Most of this gas remains locked up 
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Fig. 15.— Evolution of star formation density with respect 
to the total computational box for all galaxies that have M^c > 
10^ h~^ Mq at each redshift. This is different from Figurelm which 
provides star formation histories for g alaxies chosen at 2 : =^- Th e 
dashed line is the best fitting line from|Madau & Dickinson|(|2014|), 
and is based on UV a nd IR data. The observational points are trom 
[Bouwens et al.| ( |2015| > — blue points are uncorrected for dust, while 
red points have been corrected for the dust extinction. 

within their DM halos, cools down, and is accreted by 


the galaxies (e.g., Keres et al. 2005| Dekel & Birnboim 
Romano-Diaz et al.|2UI4) . As a consequence^ these 
galaxies become very metal-rich early-on (Fig.^, en¬ 
hance their SFRs (Fig.[^, consume their gas very quickly 
(Fig-i, and become gas-poor and very compact already 
hy z ^ 6. An overall steady decay in the gas fraction 
within the NW galaxies is correlated with their SFRs 
(Fig.§. The smallest effect is observed in the lowest 
mass-bin, i.e., 7.5 < logM^^^ < 8.5. This happens be¬ 
cause such objects are less prone to the feedback effects 
due to their respective, already very low SFRs. 

Our conclusions regarding the NW model pertain 
mostly to the decline in the gas fraction in intermediate- 
mass and massive galaxies, and consequently the decline 
in their respective SFRs. This implies reddening of such 
galaxies due to aging of their stellar populations at rela¬ 
tively high-z (e.g., 2 : > 3). However, such a trend seems 
to be in contradiction with the observed rise in the cos¬ 
mic SFR which reac hes its peak at z ^ 2 — 3 (see Madau 


& Dickinson (2014) for the latest review), although the 
latter SF'R is averaged over all possible environments. 
Furthermore, this model is in disagreement with the mea¬ 
sured SFR at z ^ 6 —8 (Fig.[^ due to an early gas deple¬ 
tion, and deviates most from the empirical sSFR - 
median shown in Figurej^ 

For a more direct comparison with observations of 
high-z galaxies, we have calculated the global star forma¬ 
tion density in our computational box for galaxies with 
> 10 ^ Mq (Fig. 15). We include all the wind 
models in the overdense and normal regions. Further¬ 
more, we have adde d the best-fitting curve from [Madau 


& Dickinson (2014), given by their Eq. 15 (see also their 
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Eig. 9), as well as av ailable observational points from 
Bouwens et al.j (2015|). 

A number oi conclusions can be made based on Eig- 
nrep^K) Eirst, the overdense models lie well above the 
normal density models. However, their slope is in good 
agreement with the observational data points, while 
their amplitude is above that of the data, as expected 
(|Romano-Dfaz et al.||2014| |Yajima et al.||2015|). Second, 
an extrapolation of tne trend from the UCW and UVW 
models to higher 2 : show that their respective curves have 
a slope which agrees well with the observational points at 
2 : ^ 7—10, although it is somewhat steeper for 6 < ^ < 7. 
Lastly, this agreement is better than the one provided by 
the best-fitting curve (dashed line) at z >8. Hence, our 
UGR curves, although lying above the dashed curve at 
z ^ 6 — 7, provide a better match to observations at 
higher redshifts. Note, that latest corrections to the ob¬ 
servational points (both dust corrected and unc orrected) 


move up ste adily, e.g., when comparing t hem in Bouwens 


et al. (2014) with [Bouwens et"aL (2015). Given the un- 


certainties of the current observations, our small compu¬ 
tational box, and the absence of any fine-tuning or cali¬ 
bration in our simulations, the match between our UGR 
models and observations is acceptable, which is impor¬ 
tant when one discusses the overdense regions. 

An immediate consequence of the NW model and the 
locked-up gas content within their respective DM halos at 
2: > 6 is related to the reionization process of the universe 
and the metal pollution/enrichment of the IGM. If the 
metal-rich gas is locked within the DM halos and is highly 
concentrated, the dust content within those halos and 
galaxies will be relatively high, making difficult for the 
UV photons to escape and ionize their surroundings. The 
escape fraction of photons will be lower than the expected 
fraction of ^ 0.3 required to maintain the universe fully- 
ionize d by ^ ^ 6 (Einkelstein et al. 2012 Kashikawa et al. 
2011), resulting in a delay in the process of reionization. 

In comparison, the GW model points to a relatively 
high temperature IGM gas (Eig.|^ that extends well be¬ 
yond the halo boundaries, deep into low density regions, 
where most of the IGM has reached T>2 — SxlO^Kby 
2 : ^ 8. This temperature can collisionally ionize the HI, 
even in the absence of UV photons. This is in contradic¬ 
tion with the currently accepted scenario for reionization 
in which the neutral gas in the IGM is photoionized by 
the UV photons coming from galaxies and quasars. By 
2 : ^ 6, the GW model shows an overheated IGM gas, and 
we observe the formation of a hot bubble constituted of 
a shock-heated gas with a temperature up to r\j 10^ K 
in the central region. This GW model gas, potentially 
detectable as an X-ray bubble in this overdense environ¬ 
ment, is absent in its UCW counterpart (Eig. 14). Impor¬ 
tantly, both GW overdense model and the UUW model 
exhibit large volume-filling ionized gas with T > 10^‘^K, 
which contradicts the re-ionization constraint. 

Both NW and GW models exhibit a minimal depen¬ 
dence on the host galaxy properties. However, galactic 
outflows are expected to be driven by various mecha¬ 
nisms, including radiation pressure, S N and stellar winds, 
and photo-heating of H H regions (e.g., Choi fc Nagamine 


2011 ][Hopkins et al.|2012||Agertz & Kravtsov|2U15|). big- 
nihcance of each process depends on the galaxy proper- 
ties and one should consider all mechanisms in the galaxy 
formation simulations to properly account for the cumu¬ 


lative effects. Within this framework, the galactic out¬ 
flows are hardly described by a simple scaling relation¬ 
ship. Consequently, we turn to the eff ects of a more 
complex type of models, the VW (section |2^ . 

In the VW model, the winds properties depend on 
the SER of the host galaxy, calibrated at \ow-z. Al¬ 
though this is a rather simplistic parametrization of a 
very complicated process, it introduces self-regulation in 
the galaxy growth. This growth is based on the local, in¬ 
stantaneous SER rather than being independent of host 
galaxy properties, as in the GW model. 

In our simulations, the VW model exhibits the most 
complicated (e.g., nonmonotonic) behavior in terms of 
interplay between galactic outflows, gas consumption and 
SE among all the wind models. This behavior is espe¬ 
cially pronounced for the gas fraction, /gas,gal- On the 
average, the gas fraction in VW galaxies lies between the 
gas-poor NW galaxies and the GW gas-rich ones, at the 
end of the simulations. This is also true for the median 
SERs. The VW objects also show a strong decline of 
/gas,gal with the galaxy mass, in a sharp disagree¬ 
ment with the GW models. 

When comparing with observational data (e.g., Eig- 
urel^, a case in favor of the VW prescription can be 
made. Due to the direct coupling between galaxy and 
wind properties imposed in the VW prescription, the 
SER is correlated with the mass loss rate. This self¬ 
regulation is very important for the intermediate mass 
galaxies, 10^'^ < < 10^^, which are the most com¬ 

mon population at these high z, and probably the main 


contributors to reionization, but see Madau & Haardt 


(2015). We find that VW galaxies exhibit a tight corre 


lation sSER - which lies within Icr from the observa- 
tionally determined median shown in Eigurej^ 

Observational results concerning the state of the IGM 
metallicity at high- 2 : make it difficult to draw conclu¬ 
sions in favor of specific wind models in identical environ¬ 
ment, based on the differences between their efficiency in 
spreading metals on large-scales. Recent studies focusing 
on the statistics of GIV absorbers in spectra of distant 
quasars have shown a decline in their occurrences fo r 
2 : > 5.3 (Becker et al. 2009 Ryan-Weber et al. [2009). 
Ryan-Weber et 2009[ ) found a cosmic mass density 

of GIV ions at 2 : ^ 5.7 of Uciv = 4.4 x 10“^ 


of 3.5 lower than the corresponding value for lower red- 
shifts, and derived a corresponding lower limit for the 
IGM metallicity of Ziqm ^ 10“^ Zq. 

The enrichment levels in the IGM at 2 : ^ 6 in all of our 
simulation runs seems to satisfy this criterion based on 
the IGM metallicity maps shown in Eig.|^ Even the NW 
model is able to rea ch peak metallicities of Z ^ 10' -3^0 
in the IGM (Eig.I^ despite the clear lack of metal spread 
on laxge-scales and a very small volume filling factor 
(Eig.jM. The reason for the relatively efficient IGM pol¬ 
lution even in the NW is likely caused by the accelerated 
evolution due to the presence of the overdensity which 
boosts SER and, consequently, the metal production in 
this region. 

It is worth noting that other simulations used different 
outfiow prescriptions to model the observed IGM metal¬ 
licity at low 2 : (e.g. Oppenheimer et al. 2QII). These 
of course do not include the no-wind niodel. [Oppen-[ 
Iheimer et al.l found that GIV absorbers are associated 
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with 10^ K gas in halos for which, in our models, this 
material is found at similar metallicity levels (Fig.[l^. 

The metal content of the IGM is affected as well by the 
presence of the overdensity because of the development 
of massive galaxies. This is in contrast with the average 
density (UCRs) models where no metal-rich outflows are 
observed. Nevertheless, we do find that the peak metal- 
licity attained in few regions in the IGM of the UGW run 
(visible as high-metallicity spots in Fig. 10) is on a sim¬ 
ilar level as the GW model. So, even though the galaxy 
population is less developed, the metallicity content is 
of the same order, albeit in limited volumes — a direct 
consequence of galaxy evolution. When averaged over 
the computational volume, the IGM metallicity in UGR 
models is lower than the GR models, with exception of 
NW. 

Heating the IGM to T > 10^ K will destroy the dust 
and facilitate reionization at much earlier times, being in 
contradiction with the latest reionization time estimates 
(jBecker et al.||2015t |McGreer et al.||2015|). Furthermore, 
the HI opacity would decrease rapidly from even ear¬ 
lier redshifts , in contradiction with the results from |Fan| 


et al. (2006). At the same time, the IGM metallicity is 
observed to be very high in our GW models, nearly sub¬ 
solar, which should leave an imprint in the spectra of the 
high- 2 : objects, such as quasars and Lyman a emitting 
galaxies. However, this has not been detected yet. Fur¬ 
thermore, it will also help the SF process at later stages, 
favoring the formation of proto-galaxies in the surround¬ 
ings in large numbers. Th erefore, galaxies, s uch as the 
so-called GR7 galaxy (e.g. Sobral et al. 2015), could be 
found in vici nity of similar objec ts, forming proto-cluster 
regions (e.g., |Trenti et al.||2012|). 

Although one could state that these are mere results 
coming from a high-density environment, this situation is 
not very different when comparing the metallicity levels 
between average and high-density environments. Their 
respective overall metallicity profiles lie within 0.5 dex, 
but their distributions for the different media (i.e., the 
IGM, halo gas and the ISM) are very similar (Fig.p!^. 
Their overall temperature profiles are similar as well 
(Fig.[^. 

Overall, the temperature in the VW model lies within 
10^ K and metallicity at \og{Z j Zq) « —2, a much more 
acceptable range which is also in agreement with reion¬ 
ization. Figurej^ shows that the VW model is able to re¬ 
produce both the observed sSFR evolution between 2 ; ^ 8 
and z 6 and the sSFR-stellar mass relation at 2 : ^ 6. 
The GW and NW models fail to satisfy simultaneously 
these two constraints. Moreover, the NW run shows a 
clear disagreement with observed levels of star formation 
atz^S — 6. Atz^G, the scaling of sSFR with stel¬ 
lar mass is underestimated ove r the entire mass ra nge 
lO^-S _ 10 ^^ h~^ Mq spanned by Salmon et al. (2015) ob¬ 
servations (Fig.l^ left). The sSFR values are also off at 
2 : ^ 8 and 2 : ^ y7 not only for the entire galaxy popula¬ 
tion, but even when considering galaxies only within the 
observed stellar mass range. 

Unfortunately, computational expenses prohibit imple¬ 
menting high-resolution studies in large-scale cosmolog¬ 
ical simulations. The full cosmological simulation still 
needs to rely on a rather simple scaling. As we discussed 


r ely c 

in section [2T^ the VW model possesses one important 
feedback parameter — the wind velocity dependence on 


the host galaxy SFR, which constrains the wind prop¬ 
erties better than the GW model. Although its value 
is tuned-up to the low- 2 ; observation ( Martin|[2QQ5 ), the 
median VW speed appears lower than that of the GW. 
As a result, the VW does not overheat the IGM com¬ 
pared with the GW, while it regulates the SF and galaxy 
growth, and shows an agreement with the sSFR obser¬ 
vations. It implies that the scaling relation that the VW 
adopts still holds in galaxy evolution at 2 ; > 6 and the 
predictions from the VW simulation, such as the galaxy 
mass function, SFR, sSFR, and /gas,gai would be more 
reliable than those in GW. 

The low SF levels at 2 : < 7 in the NW case are caused 
by the absence of winds which allows efficient consump¬ 
tion of gas by SF at high- 2 ; and gradually depletes galax¬ 
ies of their gas. This also causes metallicity in galaxies 
to reach very high values by 2 : ^ 6 for a broad range 
of stellar masses (Fig.l^. In particular, for low-mass 
galaxies, we find unrealistic metallicity levels in the stel¬ 
lar and gaseous components, causing the resulting mass- 
metallicity relation to be nearly flat in the NW case, as 
opposed to increasing with mass as shown by the other 
two wind models. Furthermore, the GW model is able to 
reproduce the observed sSFR but fails at capturing cor¬ 
rectly the decreasing trend of sSFR with stellar mass, as 
it shows a positive correlation between those two quan¬ 
tities. in addition, the GW model appears in contra¬ 
diction with observations at lower redshifts, such a s an 


excessive metal enrichment of low density gas (e.g ., Op- 
2008 Choi fc Nagaminq|2011 ). We 


penheimer & Dave 
therefore, conclude that the V W wind model more closely 
reproduces galactic wind, IGM- and halo-gas properties 
at high z, although a thorough analysis using radiative 
transfer is needed to confirm this conclusion. 

In summary, we find that galaxy evolution at ^ 6 — 12 
depends strongly on the type of feedback mechanisms 
applied. While specific prescriptions for galactic winds 
allow us to reproduce some of the observed properties 
of galaxies at these redshifts, the best fit is attained 
by using the VW model, which couples wind properties 
to the ongoing star formation and to the depth of the 
galactic potential well. No perfect wind model can be 
decided at this moment, perhaps a new wind model is 
necessary which is more tightly related to the underlying 
galaxy properties. The VW model, although quite sim¬ 
plistic, appears best equipped to reproduce the galaxy 
properties, both in overdense and average regions in the 
universe. Future improvements of galactic wind models 
must include dependency on morphological and radiative 
properties of host galaxies. 
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